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Abstract 



Contrary to popular belief, the Zeeman effect can be treated exactly in single-electron systems, for 
arbitrary magnetic field strengths, as long as the term quadratic in the magnetic field can be ignored. 
These formulas were actually derived already around 1927 by Darwin, using the classical picture of 
angular momentum, and presented in their proper quantum-mechanical form in 1933 by Bethe, although 
Y^ ' without any proof. The expressions have since been more or less lost from the literature; instead, 

,—( , the conventional treatment nowadays is to present only the approximations for weak and strong fields, 

^ ' respectively. However, in fusion research and other plasma physics applications, the magnetic fields 

applied to control the shape and position of the plasma span the entire region from weak to strong fields, 
and there is a need for a unified treatment. In this paper we present the detailed quantum- mechanical 
derivation of the exact eigenenergies and eigenstates of hydrogen- like atoms and ions in a static magnetic 
field. Notably, these formulas are not much more complicated than the better-known approximations. 
Moreover, the derivation allows the value of the electron spin gyromagnetic ratio gs to be different from 
2. For completeness, we then review the details of dipole transitions between two hydrogenic levels, and 
calculate the corresponding Zeeman spectrum. The various approximations made in the derivation are 
also discussed in details. 
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1 Introduction 

At the end of the 19th century, when Maxwell's theory of electromagnetism had been established, it was 
known that electromagnetic radiation is produced by oscillating electric charges. The charges that produce 
light were however still a mystery. At the same time, it was generally believed that an electric current was 
made up of charged particles. Lorentz suggested that atoms might also consist of charged particles, and 
that the oscillations of these particles inside the atoms might be the source of light. If this were true, then 
a magnetic field ought to have an effect on the oscillations and therefore on the wavelength of the light thus 
produced. In 1896 Zeeman, a pupil of Lorentz, demonstrated this phenomenon 1], known as the Zeeman 
effect, and in 1902 they were awarded the Nobel Prize. 

Lorentz was thus one of the first to predict the existence of the electron, which, within a year later, 
was discovered by Thomson, at least as a free particle. This, in combination with the results obtained by 
Planck and Einstein regarding the nature of black-body radiation and the photo-electric effect, led to the 
subsequent emergence of Bohr's atomic theory. 

A full understanding of the Zeeman effect can only be obtained from quantum mechanics, but Lorentz 
was nevertheless soon after the initial discovery by Zeeman able to produce a simple classical theory which 
predicted certain aspects of the polarization of the different spectral components. These were later verified 
experimentally by Zeeman. Soon after the discovery of the so-called anomalous Zeeman effect (the weak- 
field limit), around 1905, came the development of the Lande vector model of the atom. Using this model, 
with its semi-empirical rules, basically all aspects of the observed Zeeman spectra could be explained and 
predicted '2|. This was however most fortuitous, since the classical Lande g-factor happens to be identical 
to its quantum-mechanical counterpart. 

Based on the Lande vector model of angular momentum, in combination with the wave-mechanics of 
Schrodinger, Darwin 3] was in 1927 able to solve "the problem of a spinning electrified body moving in a 
central orbit in a magnetic field [. . . ] by the method of the wave mechanics in spherical harmonics" . The 
formalism presented there bears very little resemblance to how the problem is formulated nowadays, but the 
treatment is actually equivalent to an exact solution of the Schrodinger equation for the Zeeman effect. 

It is certainly amusing to read Darwin's statement that "The proper attack on this problem would 
undoubtedly be by way of the recent work of Heisenberg [. . . ] but this theory is still in the making, so it has 
not been practicable to apply it here" ; the article goes on to say that the work of Heisenberg " could give 
all the results of this paper; but it would have been harder to follow because the matrix methods are not so 
easy for most readers as are spherical harmonics." 

Today the situation is rather the opposite, and we shall therefore in this paper present the exact eigenval- 
ues and eigenvectors for the Zeeman effect in single-electron ions (i.e. hydrogenic systems), using the matrix 
formulation of quantum mechanics. These expressions are not new; they appeared as early as 1933, due 
to Bethe |j|, and they have also been published in the well-known book by Bethe and Salpeter 5], where 
several aspects of the fundamental properties of the Zeeman effect, as well as certain approximations made 
in the treatment, are discussed in depth. 

However, in neither of the above quoted works is given any indication - or reference - of how to derive 
these exact formulas. Due to this, perhaps, the presentation of the Zeeman effect found in modern text-books 
on atomic physics (see e.g. Ref.|a, Section 5.2), is limited to the discussion of the weak and strong field limits, 
known respectively as the anomalous and normal Zeeman effect, and perhaps some particular intermediate 
cases, where the Zeeman effect is also known as the Paschen-Back effect. (By strong and weak fields, we 
mean that the energy contribution from the magnetic field is large or small compared to the spin-orbit 
interaction.) 

For many-electron atoms, this is naturally the proper attack on the problem, since an exact solution is 
unattainable. However, for the particular case of hydrogenic atoms and ions, the exact treatment is not 
much more complicated; in fact, the resulting expressions for the eigenvalues and eigenvectors are quite 
simple. The spectral properties of single-electron ions have recently received rekindled interest due to a 
variety of applications within plasma spectroscopy (for a review, see Ref. 13 and references therein). In e.g. 
fusion research, strong magnetic fields are used to control the shape and position of the plasma. A detailed 
understanding of the influence of these fields on the spectral lines is crucial for a correct interpretation of 
the observed line-widths, which are employed for temperature measurements j^, |9|, among other things. 
Moreover, the magnetic field strengths employed often span the entire scale from weak to strong fields, and 
there is therefore a clear need for a unified treatment of the Zeeman effect, valid for all fields. 

The structure of this paper is as follows: in Section[21we will show how to solve the Schrodinger equation 
exactly for a hydrogenic system in a magnetic field, under the assumption that the quadratic term in the 



magnetic field strength can be ignored. We then show in Section O how expansions of the exact expressions 
reproduce the traditional results of weak and strong fields. As a bonus, a quantitative expression is found for 
determining the validity of these approximations. States with zero angular momentum must be considered 
separately, and actually correspond to an extreme case of the strong-field approximation, as discussed in 
Section^] Section|3contains a review of the dipole transition matrix elements and selection rules, and finally 
we discuss in Sectional the approximations made during the derivations. 

2 Exact solution of the Schrodinger equation 

As pointed out above, the traditional treatment of the Zeeman effect in most text-books is limited to the cases 
where the influence of the magnetic field can be separated from the effects of the spin-orbit interaction. In the 
weak-field limit, the magnetic interaction is introduced as a perturbation to the spin-orbit eigenstates, and 
vice-versa in the strong-field case. In this picture, all perturbations are already diagonal and the formalism 
is trivial. 

In the general case, to be presented here, these two effects must however be treated simultaneously. Our 
method will be to diagonalize the combined perturbation in a basis set consisting of the solutions of the 
field-free, non-relativistic Schrodinger equation for a single-electron ion. Depending on the orbital quantum 
number i of the states considered, the matrices will be of various dimensions. A significant simplification 
will however emerge from the algebra, and the final energy matrix will separate into a number of 2 x 2 
sub-matrices which are easily diagonalized. The case ^ = 1 is treated in some detail in Ref. l(j, from where 
the idea for the method has been adopted and further developed. 

2.1 The unperturbed system 

The non-relativistic Schrodinger equation for a hydrogenic system in vacuum, ignoring for now the spin-orbit 
interaction, is 

^ -' ^^' ^ ^W = H^^r) = Eot^ir), (1) 
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where fj, = mcM/{mc + M) is the reduced electron mass (me and M are the electron and nuclear masses, 
respectively) and Z the nuclear charge. Because of the spherical symmetry, the wavefunction ip can be 
written as a product 

\n£mi) = tlJntmt{r,0,ip) = Rnt{r)Yim{0 , ^) (2) 

of an angular part, given by the spherical harmonics YimiO, (f), and a radial part Rni{r). 

The energy of each eigenstate | n£me) will depend only on the principal quantum number n as 

where a — e^ /Aireohc « 1/137.036 is the fine-structure constant. Introducing rclativistic corrections for the 
momentum operator adds another term 



A x^ 1 9 (Za) 



2 " n 



4 



3 



4 ^-fl/2 



(4) 



to the energy, breaking the degeneracy in the i quantum number. The degeneracy in the m quantum number 
is only broken through an external field such as the magnetic field about to be studied. 

It is not entirely obvious which value of the electron mass to use in Eq. |@J) , but since the perturbation 
in itself is a small correction, it will make an insignificant difference whether we use the electron mass rUc 
or the reduced mass /i |6|, footnote on p. 196]. Therefore, we shall in what follows always use the standard 
electron mass in perturbation theory. 

2.2 Interaction with an external field 

The interaction between an atom and an external electromagnetic field is introduced through the vector 
potential A. From Maxwell's equation V • B = (there are, probably, no magnetic monopoles), it is clear 
that one can choose to represent any magnetic field B as B = V x A. The choice of A is not unique, but one 



can generally write A = (B x r)/2 + V^, where ^ is an arbitrary scalar potential. Note that over atomic 
distances, the macroscopic magnetic field can be assumed to be constant in magnitude and direction. 

The non-relativistic Schrodinger equation (relativistic corrections and other approximations are discussed 
in Section ISJ for a spin-less electron in an external magnetic field can be derived from the Klein-Gordon 
equation 5], and is 
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where we also have used the freedom of gauge invariance to set V • A = (the Coulomb gauge) , which means 
that the scalar field ^ must be a constant, which we arbitrarily may choose as zero. 
The term linear in A can be written 

A V = :^B L, 

rric Zrric 

where L — —ih{r x V) is the angular momentum operator. Similarly, the quadratic term becomes 
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Comparing the magnitudes of the linear and quadratic terms, one finds jg] that their ratio is of the order 
B ■ 10^^/Tesla for atomic systems. Except for extreme cases, such as neutron stars, the quadratic term is 
hence completely negligible, which simplifies the mathematics a great deal. (When included, the quadratic 
term will lead to diamagnetism |5|.) For very strong magnetic fields, when the spin-orbit interaction can be 
neglected altogether but the quadratic term not, the energy levels can be found by different methods |ll| . 

So far we have used a semi-classical picture of the electron by ignoring the spin, which however contributes 
an additional, intrinsic, angular momentum. To the magnetic moment of the electron from the classical 
angular momentum 

ML = --^L^-^iBL/h, (6) 

ZrUc 

where fiB is the Bohr magneton, we should therefore add a contribution from the spin 

Ms^~gsT^S^-gsfiBS/n. (7) 

According to the original Dirac theory, the magnetic moment of the electron due to its intrinsic spin is 
determined by gs = 2. However, using quantum electrodynamics (QED), one finds that the value of the spin 
gyromagnetic ratio gs is not exactly integer, but 

Os a 0.328a2 , , 

2±^l + — + ^ + ...« 1.0011595. 8 

This correction is called the anomalous magnetic moment of the electron. Since it poses no problems to 
retain a general value of gs in our calculations, we will do so. Moreover, this allows the formalism to be 
extended (with proper modifications, accounting for the relative dielectric constant and the electron effective 
mass) to e^. solids, where the values of the gyromagnetic ratio may differ significantly from 2, and even be 
negative |l^ . 

The energy of a magnetic dipole - which the electron effectively becomes due to its magnetic moment - 
in a magnetic field is 

i/;,g„ = -(Ml + Ms) • B = ^ (L + 5,8) • B, (9) 

and we now also introduce the energy from the spin-orbit interaction 

H',, = ^r) (L . S) , (10) 

where for the one-electron central-field case 

Zm^c' 47reo r'^ 

In the field- free Schrodinger equation, Eq. |^, the geometry is spherically symmetric, whereas the mag- 
netic field defines a particular direction of symmetry. We make this direction the z-axis of our system, which 
reduces the scalar product in Eq. to B • (L + g^S) = B [L^ + gsS^), where B = |B|. 



The total non-relativistic Hamiltonian for the electron in a hydrogen-like atom in an external magnetic 
field thus becomes 



H - Ho + H' + H' 



magn 



We assume that the relativistic correction Ai^ioi is still given by Eq. I^J in the presence of the external field. 

2.3 Perturbed energies and eigenstates 

To diagonalize the perturbation 

H' = H'so + ^magn (13) 

we first need to introduce a proper set of basis states. This is actually not trivial |5!, Section 45/3], but to 
simplify matters we shall use the Pauli approximation, in which each state is a two-component spinor and 
the spin operators are described by the Pauli spin matrices. 

We can therefore use as basis the eigenstates of the unperturbed Hamiltonian Hq, arranged in a proper 
spinor notation. In the unperturbed system, the good quantum numbers are, in addition to n and £, also 
mi and tus, although the eigenstates are degenerate in the m:s. In the presence of the external field, the 
individual m:s are no longer good quantum numbers, but only their sum m = me + m,s. 

We can therefore anticipate that each new eigenstates will be a linear combination of exactly two un- 
perturbed eigenstates, corresponding to nis = i^ and m,e chosen accordingly But since the spin-orbit 
interaction is not diagonal in the unperturbed eigenstates, we need to carry out the derivation using the full 
expansion in the unperturbed eigenstates 

\n£m) — 2_^2_^ Inlmim^s) {n£mim,s\n£m) , (14) 

where TOs and m^ take on all allowed values. To avoid unnecessary notation, the indices n and £ will be 
suppressed most of the time, as all calculations are performed for fixed values of these quantum numbers. 
In all matrices below, the basis states are assumed to be arranged in the following order: 

(m,,m,) = (^,t), (^,i), (^-1,T), (^-l,i), ..., (-^,T), (-^,i) 

where up and down arrows mean rn^ = ±^ respectively. 
The contribution to the energy due to H' is 

AEnem(B) = {n£m\H'\n£m) 

= W E ^n£m\mem.) Um. ^ + ^1^/3 m',m'\ ^m'.m'Mm) , ^''^ 
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where we have defined 

l3 = B/Bo, (16) 

with Bq — W//iB where 

The expectation value W is not defined for ^ = 0, since then there is no spin-orbit interaction, but instead 
an energy contribution named the Darwin term. We will consider that case separately in Section^ and may 
assume £ 7^ in what follows. One can then show 6] that for the hydrogenic eigenfunctions. 
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The magnetic field term, 
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is already diagonal in the present basis, but the spin-orbit interaction Jj ■ S — L^Sx + LySy + LzSz is not. 
Since the orbital angular momentum operator L and the spin angular moment operator S act on separate 
vector spaces, spin and position space respectively, they trivially commute. Therefore, we can separate space 
and spin, so that 
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This expression corresponds to an outer or Kronccker product between the space and spin matrices. This is 
defined as 
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(22) 

from which generalizations to larger matrices should be obvious. Each element in the product matrix is a 
product of one element in the first and one in the second matrix, spanning all possible combinations. Within 
each marked sub-matrix, one finds the matrix B multiplied with the subsequent elements in A. 
The required matrix elements for the spin part are immediately given by the Pauli spin matrices 
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whereas for the space part we make use of the ladder operators 



{L+ + L^)/2, 
{L+-L^)/2i. 



(24) 



whose action on the unperturbed eigenstates is defined through the relations (with conventional phase 
choices) 

L± \e, mi) = h^iit + 1) - mi(mi ± 1) \i, m^ ± 1) . (25) 

To write down the matrices corresponding to the operators Lx and Ly is somewhat tricky, since the 
sizes of the matrices depend on the value of i\ the dimensions will be 2i? + 1, i.e. the number of possible 
m-values. The 2-component is however diagonal, since the operators Lz and Sz correspond to the good 
quantum numbers rat and nis of the basis eigenstates, so we can write it generally as 



{mi\Lz\m!^ {ms\Sz\"m!^ = mimgh Smim.\^m,m'^- 



(26) 



This will be the only contribution to the diagonal of the spin-orbit term, since the ladder operators L± (and 
hence Lx and Ly) do not couple states with the same mi or m^. 

Moreover, although each of the matrices {mi\Lx\m,g) and {me\Ly\m'i) contain several off-diagonal ele- 
ments, an important simplification takes place when they are added together to form (m^ms |L • S|?tt-^toJ,), 
and the matrix separates into a sequence of 2 x 2 sub-matrices, plus the corner elements which are uncoupled. 
Clearly, this represents what was predicted earlier: only states of equal m = mi + rn^ are coupled. 



Using these two observations, it becomes possible to present the structure of the spin-orbit matrix as 



TTl^TOs 



2L-S 



(. 





: 







: 





-(. 


? \ 







: 





7 


l-l : 







: 








: 


-{i- 


1) 


? \ 








: 


? 




i-2 : 















(27) 



where the so-far undetermined elements are denoted "?". The structure is emphasized by the dotted lines, 
which enclose 2x2 areas in which all matrix elements correspond to given quantum numbers m and m'. It 
should be apparent that the non-zero matrix elements will be confined to the block-diagonal, and that in 
each such block the relation m — m! will hold. 

The corner elements of the matrix are trivial, so we focus on the interior of the matrix. To find the 
unknown matrix elements, we may use the same matrix but with a different division into sub-blocks: 
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This division corresponds exactly to the direct matrix product H22II . and we can therefore easily trace the 
contributions to the unknown elements. 

First consider the top right question mark (boxed) in Eq. \2T\ . In this area of the matrix, ra'g = £ — 1 
and mg — £. These states are coupled by L-|_ and we get for the sub-matrix 



r /o 
2 I? 



LyGy 



{me ^ l\L+\m'p = £ - 1) f -a^ + —cfy 
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For the unknown element in the third row, second column, the only difference is that the involved states 
are coupled by L_ instead. The prime is only a book-keeping label; mi and m^ are interchangeable as long 
as i_|_ is replaced by L_ and vice versa, and hence this second unknown clement is exactly the same as the 
first one, with the only change of removing the primes. 

In the positions of the matrix where the question marks appear, we always have nis = 1/2 or m'^ = 1/2. 
Therefore, on these positions, mi ^ m — 1/2 (with or without primes). These observations are however not 
particular for this sub-block of the matrix, but apply generally, and thus all non-zero off-diagonal elements 
can be written on the same identical form y^£{£ -I- 1) — (m, — 1/2)(to -|- 1/2), since furthermore m, — m,' 
within each sub-block on the diagonal, as noted above. 

Finally we rewrite the diagonal elements, which have contributions from both the LzSz matrix elements, 
Eq. if^ . and the magnetic interaction, Eq. it^ : 
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When rus = ±1/2 we must have mi = m^ 1/2 and thus 

V{mi, = TO =F 1/2, TO^ ^ ±1/2) = ± ( m =F t: 
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Again, these expressions will be valid anywhere on the diagonal, and we have thus proven that all the 2x2 
blocks along the block diagonal can be written on a common form. 



-{m+l/2)+(3[m-{gs-l)/2\ ^ i{i + 1) - (m - l/2)(m + 1/2)' 



v/£(^ + 1) - (m - 1/2) (m + 1/2) m - 1/2 + /3[m + (.g, - l)/2] 



(29) 



valid for all allowed values of vn, except of course the corner elements \m\ = £ + 1/2. These are however 
trivial, and give the eigenvalues (cf. Eq. H15|) ') 
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To find the other solutions, we need only diagonalize the matrix H29() . which has eigenvalues 

W 
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2m/3 - 1 ± V/?^(3s - 1)2 + 4m/3(g^ - 1) + {21 + 1)^ 
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which thus is the total energy shift due to the spin-orbit interaction and the magnetic field. The correspond- 
ing eigenstates can be expressed in terms of the original \mi, mg) states as 



\m)^ = sinif \m+ 1/2, i) + cosip \m — 1/2, j) , 
\m) _ = con if I'm + 1/2, 4) — simp |77i — 1/2, |) , 



where the angle ip is given by 
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The absolute signs appear here since an explicit sign convention has been introduced in Eq. (|32ab ). 
Examples of the calculated Zeeman levels based on these results are shown in Figs. Q and El 



3 Strong and weak field expansions 

The formalism as presented above makes no assumptions of the relative strength between the spin-orbit 
interaction and the magnetic field. In the traditional derivation, however, the one of the two which is stronger 
is first treated as a perturbation to the central-field Hamiltonian, which gives a set of new eigenfunctions. 
Subsequently, the other effect is treated as a perturbation to this new system, using its eigenfunctions as a 
basis set. 

In many practical cases the magnetic fields in question are either "weak" - laboratory field strengths are 
typically of the order 10-100 mT - or "strong" (superconducting electromagnets or astrophysical objects), 
making the traditional approach seemingly more appealing, since the algebra and the resulting expressions 
are somewhat simpler (but only very little!). 

It is therefore desirable to attempt to derive these limits from the general case, and also to obtain 
a quantitative criterion for when these approximations are valid, which is generally brushed over in the 
traditional texts, where "strong" or "weak" fields are merely assumed. For the weak- field case, it is necessary 
to assume gs = 2, and therefore we will also assume this for strong fields, in order to facilitate the comparison 
with expressions in the literature. 

As our starting point we will use the expression for A_E,^, Eq. H31|l . The parameter determining the 
relative strength of the magnetic interaction to the spin-orbit interaction is f3. Assuming that this quantity 
is large, i.e. (3 ^ {2i 4-1), we can neglect both {2£ + 1)^ and the term linear in /? in the square-root, and 
also the term (—1) outside, giving 

A^± = B^iB {m ± 1/2) . (34) 

The solutions to the corner elements given by Eq. H30|l also reduce to this expression if we ignore £ compared 
to (i + l)/3. This result is trivially obtained from Eq. H15|l if the spin-orbit term L • S is ignored, and thus 
corresponds to the well-known expression for the normal Zeeman effect. 

In this strong- field limit, the states are split in m but degenerate in £ (except of course for the dependence 
on £ in the relativistic correction, Eq. Q), as seen in Fig.|21 The eigenstates are particularly simple, since the 



denominator in Eq. (|33|l becomes very large, and so ip ~ 7r/4. Since the levels tend to group, and because of 
the selection rule Am = 0, ±1 to be discussed shortly, the spectrum for normal Zeeman effect of a transition 
n — > n' will consist of three equidistant lines - a Lorentz triplet 6, p. 211]. 
In the opposite limit, where /3 ^ {2i +1), we get 

1 4to/3 




^^m ~ — 2m/3 - 1 ± (2£ + 1) <j 1 
Wi + BfiBm 



2(2£+l)2 
2^ + 2 (35) 
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2£ 
-Wie+l) + BfiBm^^-^, 

if we neglect the term quadratic in B. The solutions to the corner elements are also incorporated in this 
formula (the top expression), which is easily shown - without approximations - from Eq. (j^OJ), if we realize 
that m/{£+l/2) = ±1. 

We can identify the first term in Eq. (|35|l as the contribution from the spin-orbit interaction if the 
magnetic field is ignored. The two cases correspond to j = £ ± 1/2, where j is the quantum number 
associated with the total angular momentum J = L + S which is conserved in the absence of a magnetic 
field. The energy contribution from the magnetic field can now be evaluated, within this approximation, 
through the expression Q 



T«^' + ^*)^T 



Ai^magnK) = ^^(i. + 25, ) - ^B(J, 



(this is where we need gs = 2). The expectation value of J, is hnij by definition, and one may show that 

{nisjmj\Sz\n£sjmj) — {g — l)m,jh (36) 

where s is the spin (i.e. 1/2 for an electron) and 

j{j + l) + s{s + l)-£{i+l) 
2jU + 1) 

is the Lande factor. Thus A_Eniagn = gfJ'BBmj, which can be shown to give exactly the remaining part of 
Eq. (|35(l . We now also realize why the corner solutions ended up in the top expression of Eq. 135(1 . since 
\ni\ = £ + 1/2 only is allowed for the triplet state j = £ + 1/2. 

In the weak-field case the distance between two levels corresponding to subsequent values of m within 
the same multiplet will be the same for all levels (cf. Fig. ^. This case is, for historical reasons, called 
anomalous Zeeman effect, but is actually the more commonly encountered one. Since the splitting of the 
levels for different multiplets is different (i.e. dependent on j), there will be more lines in the spectrum than 
for the normal Zeeman effect. 

The approximate expressions for strong and weak fields can thus be derived as expansions of the exact 
result, Eq. H31() . We may now also specify quantitatively under what conditions these approximations hold, 
since this is determined by whether /3 — B/Bq = Bfis/W is much greater or much smaller than 2£ + 1. 
Inserting the constants into Eq. 1)17(1 gives 

yi yi 

W = 3.62261 • lO"'^ eV = 2.92183 cm"\ 

7 7 

Ba = 6.2584 Tesla 
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where 7 = ti^£{£ + 1/2) (^ + !)• Thus the approximations are valid if B is much larger or smaller than 

(Za)4 niec^ Z^ 



nH{£+l) 2^B n3£(£+l) 



12.5 Tesla. (38) 



For light atoms and states with n and £ less than 4, this expression takes values between 0.1 and 10 Tesla; 
for heavier atoms the fourth power of Z will dominate and, except for high shells and strong fields, the 
weak-field approximation will almost always be valid. 



As an example, let us consider the states n = 3 and 4 in HcII {Z = 2). Inserting £ — 1 and 2 for n = 3 
(for i — the strong- field approximation is exact; see the following section) the threshold field strengths 
come out as 3.70 and 1.23 Tesla respectively. For n=4 the corresponding values for i'=l, 2 and 3 become 
1.56, 0.52 and 0.26 Tesla. As briefly mentioned in the Introduction, a modern application of Zeeman effect 
is to use fusion plasma spectroscopy in order to determine the plasma temperature, by e.g. studying the 
transition between precisely these states in Hell [lj|- Typical magnetic fields in a tokamak the size of JET's 
lie in the range 1-3 Tesla, and we thus clearly see that such an investigation cannot be carried out using 
either of the traditional approximations. 

4 Treatment of s-states 

As mentioned above, the expectation value W is only defined for £ ^ (cf. Eq. (|19|l 'l. But, in light of the 
previous section, we now realize that i — corresponds to an extreme case of Eq. H34() . where the magnetic 
field is infinitely much stronger than the spin-orbit interaction. Hence, without approximation, the energy 
eigenvalues are (here we may keep a general value of gs) 



Ai?± == B^IB [m ± ^^^j = ±BfiB9s/2, (i - 0), (39) 

since for ^ = we have m = nis = ±1/2. To this we should also add the Darwin term [g 

^^- = 2;^4^ (-OOl^(.)hOO) ^ ^i-^. (40) 

5 Transitions and matrix elements 

To calculate the Zeeman pattern of the spectral line from a transition between two different nil states, we need 
to consider the selection rules and the transition matrix elements. Since the eigenstates (|32ab ') of the total 
Hamiltonian, Eq. (|12|) . are simply linear combinations of hydrogenic states \nimims) ^ it is however enough 
to employ the well-known results for transitions between two hydrogenic states (derived in the absence of 
any magnetic field). 

The derivation of the matrix elements for spontaneous emission is a procedure found in any textbook on 
the theory of atomic spectra (see e.g. Ref. |a, Chapter 4). We therefore merely state the results in order that 
this paper can be used as a self-contained reference on how to calculate the Zeeman spectrum. We will then 
give an illustrative example of the resulting spectrum when these expressions are applied to the Zeeman 
problem at hand. 

5.1 Spontaneous emission matrix elements 

In the dipole approximation, the transition probability for spontaneous emission from a state 6 to a lower 
state a is found from Fermi's Golden Rule to be 

w:,{n) dn = ciiOba) \e ■ vbal'' dn, (4i) 

where dfl is the solid angle element in which the radiation is observed, and e is a unit vector in the polarization 
direction of the radiation. The overlap matrix element 

Tfja — {n' i' m'i/m's\v\nllmims) — {n'£'mi\r\nimi) Sm'^m^, (42) 

where we introduced the prime notation for the higher state b for a cleaner notation, and used the fact that 
the operator r does not act on the spin space. Hence we obtain the first selection rule, ArUs = 0. 
The pre-factor in Eq. (|41|l is 

f2zY {AEf 1 

where AE = TiLo^a = Eb — Ea- It is generally only meaningful to compare the relative intensities of Zeeman 
components for the same chemical element (same Z and /i), and hence the only factor we really need to 

10 



retain in C is the energy difference AE. When calculating the spectral line intensity from the transition 
probability, an extra factor AE enters, making the observed intensity proportional to (AE)'^. 
Introducing the spherical tensor components of e 



1 

7^' 



eo = e^ £±i'=^T—^{^x±iey), (44) 



and correspondingly for rba, we can write 

1 1 



where 



q=-l q=-l 



I\„ , = a/ — K7 f Yf*,Yi„Yem,dfl. (46) 



with the radial overlap integral 



T^if ^ I Rn't{r)Rnt{r)r' dr. (47) 



Here Rni{r) is the radial part of the hydrogenic wavefunction, to which we will shortly return. 

U sing either the Wigner-Eckart theorem or the properties of the spherical harmonics, one can show 
that |l4L Section 11-4 and Chapter 14] 



^1',^'ni^, = {-If >-'''' ^/i;'Kf r ^ M , (48) 



\rai q —irigj 

where ^> = max(^',£). The properties of the Wigner Sj-symbol now enable us to identify further selection 
rules, since the 3j-symbol vanishes unless £' = £+! and mi + q = m'^. Thus we obtain the £ selection rule 
which requires change of parity; dipole transitions only connect states of opposite parity, and the parity of a 
hydrogenic state is (—1)^. Moreover, since q = 0, ±1 we may write the other selection rule as Ami = 0, ±1. 
In summary, the dipole selection rules for hydrogenic states are 

Ams = 0, A^=l, Am^ = 0,±l. (49) 

5.2 The radial integral 

Including proper normalization, the radial part of the hydrogenic wavefunction becomes 



where >C^J^_j^(x) are the associated Laguerre polynomials, a^ = 47reo?i //^e^ is the Bohr radius (adjusted 
for finite nuclear mass), and we have also introduced p — 2Zr/a^ as the dimensionless length scale. 
Inserting this into the integral H47|) , one arrives at 

</ = ^T-J (51) 

with the dimensionless form factor 



"^ 2Y (n + £)! (n'+f)! n^+2 „'^'+2 

/•oo 

X / e-f (^+^) /+^'+^ L^:^l,{pln) eS^l^ApH) dp. 
Jo 
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(52) 



5.3 Polarization 

We now return to Eqs. H41|) and H45|) to study the polarization of the emitted radiation. Any arbitrary 
polarization can be expressed as a linear combination of two independent directions ei and 62, such as left 
and right circularly polarized, or two perpendicular linear polarizations. For this discussion we choose the 
latter. We also introduce a wavevector k of unit length (as are the e-vectors) pointing in the direction of 
propagation of the emitted light, so that the three vectors ei, 62 and k form a right-handed system (see 
Fig. 131. As before we let the z-axis be the direction of the magnetic field. Taking, arbitrarily, 62 to lie in 
the xy-plane and ei pointing downwards, we have 

k — (sin 6 cos 0, sin 9 sin cj), cos 9), (53a) 

ei = (cos 9 cos 0, cos 9 sin 0, — sin 9) , (53b) 

62 = (— sin0, cos(/>, 0). (53c) 

Next we use the properties of the 3j-symbol in Eq. H48|l . which is non-zero if and only if m'^ ~ m( + q, 
as already stated. Since each of the three values of q—0 or ±1 only may be satisfied one at a time, we have 
three separate cases to study. These are for historical reasons referred to as n and ct^ (from the German 
words for parallel and perpendicular, " senkrecht" ) . 

TV In this case Am£ — 0, and consequently wc only need to consider 1° (we suppress the quantum numbers 
in this part). This gives the transition rate 

iy,^fc(7r)=C(c^bJsin2 0|xO|' (54) 

for a photon with polarization ei, whereas the rate is zero for polarization 62. Since the k- vector 
is defined as the direction in which the radiation is observed, we find that if the emission is viewed 
longitudinally (k parallel to the magnetic field), 6* = and so the tt component is absent. In the 
transverse direction (9 = t^/'2-), i.e. in a direction perpendicular to the magnetic field, we have ei = 
(0, 0, —1) and thus the tt component will be plane polarized with e = ei in the direction of the negative 
z-axis. 

a~^ This case corresponds to Am^ = m'^ — mi = — 1. Now the only contribution will be from T~^ , and 
since the two polarization directions are orthogonal and independent, we may sum the probabilities 
instead of the amplitudes. In result, we obtain 

M/,\(a+) = C(c.,,)i±^|I-l|^ (55) 

The a^ component is plane polarized with e — 62 in transverse observation, whereas in longitudinal 
observation it is circularly (left-hand) polarized. The intensity of each a component is always half of 
the TT component in transverse observation. 

(T~ For the case Atti^ = 1 we get the analogous expression 

I^:,(a-)=c(c.,.)i±|^|I+l|^ (56) 

The <T~ component is also plane polarized in transverse observation, and (right-hand) circularly polar- 
ized in longitudinal. 

From here on it is a simple matter to calculate the intensities of all transitions in the Zeeman pattern; 
see Fig. 0] for an example. The number of components increases rapidly with the quantum numbers; for a 
7 — 8 transition the total number is 512, although many will be very weak. 

6 Discussion of approximations 

In this section we summarize the approximations made in the treatment. They are not many, as we have 
strived to keep things as general and exact as possible. For an even more detailed discussion, see Ref. |^ 
Sections 47a-e. 
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• The quadratic term 

In Scction f2.2l we ignored the term in the Hamiltonian quadratic in A, and we also showed that for fields 
of normal strengths (less than 1000 Tesla or so), this approximation is certainly valid, and simplifies 
the algebra considerably. 

A particular effect which enters if the quadratic term is retained is the possibility of mixing states with 
different i quantum numbers. The linear term in A commutes with L^, making i a good quantum 
number, but the quadratic term does not. Specifically, it introduces mixing between states with (.- 
values differing by ±2. Again, this will only be of importance for very strong magnetic fields and can 
often safely be ignored. 

• Relativistic corrections 

Using the Pauli approximation instead of the Dirac theory means we are only calculating the Zeeman 
effect to lowest order in a^ . Assuming that all perturbations (magnetic field and spin-orbit interaction) 
are small, we remedy this by using the same correction Eq. Q) as for the field-free case to correct for 
the non-relativistic momentum operator. Effects beyond this should be negligible. 

• The anomalous magnetic moment 

In the approximate treatment of the anomalous Zeeman effect (weak fields) it is necessary to use the 
value (7s = 2. As it turned out, we were not restricted by this approximation in our exact derivations 
(note the factors 5s — 1 and 5s/2 appearing in the final expressions, Eqs. %M)\ and 1)31(1 '). although the 
influence is very small - unless of course gs differs significantly from 2. 

• Nuclear motion and hyperfine structure 

We have already taken into account the nuclear motion by using the reduced mass /i of the electron 
for the central-field Hamiltonian. In the perturbations, we have kept the electron mass since the 
perturbations themselves are small, so any corrections to them are even smaller [^. 

Hyperfine structure due to the nuclear magnetic moment is in general orders of magnitude smaller than 
the fine-structure splitting (spin-orbit interaction). Therefore, we expect that hyperfine structure will 
not perturb the Zeeman spectrum except for extremely weak magnetic fields, where the Zeeman pattern 
can be heavily influenced. 

7 Summary 

In conclusion, this paper has demonstrated how the problem of a single-electron atom or ion in a static 
magnetic field can be solved exactly, under the only simplifying approximation that the term which is 
quadratic in the field strength can be neglected in the Hamiltonian. No assumption was made regarding the 
value of the electron spin gyromagnetic ratio, which allows the obtained results to be readily generalized to 
applications in e.g. solids. 

By using the results presented here, a unified treatment of the Zeeman effect becomes possible over the 
entire range of fields presently employed in e.g. fusion plasma, where the influence of the Zeeman effect on the 
plasma temperature measurements has been demonstrated to be significant in many cases. The well-known 
formulas for the weak and strong field limits are easily obtained as expansions of the general results, and in 
addition it was possible to obtain a quantitative expression for the fields at which these approximations are 
valid. 

The obtained general expressions for the eigenstates and their corresponding energies are not very com- 
plicated. As a final remark, one may therefore point out that the derivation, as presented here, has certain 
pedagogical merits and provides a good example of how to apply the power of matrix algebra to quantum- 
mechanical problems. 
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Figure 1: The total energy shift due to both spin-orbit interaction and an apphed magnetic field for the 
n — 2 shell in Hell. The energy zero-point is chosen at 2s ^Si/2 at i? = 0. The states are marked by their 
quantum number m with a following ± corresponding to the notation in Eqs. (|31|1 and H32ab ). except for 
the states with |77i| = £ + 1/2 for which there only are two unique solutions, Eq. (|30l) . For small magnetic 
fields, the states are split linearly according to Eq. (|^ . Note that without any magnetic field, the two levels 
2s ^Si/2 and 2p ^Pi/2 are degenerate (the small difference due to the Lamb shift is ignored). Also see Fig.|21 
where the field scale is extended. 
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Figure 2: In this figure, the same system as in Fig. QJis shown, but for stronger fields. As the magnetic field 
approaches high values, the pattern converges into five states, according to Eq. (|34(l . Since m takes on the 
values ±1/2, ±3/2, there are only five possible values of m ± 1/2, namely 0, ±1, ±2, and the degeneracy will 
obviously be one for the ±2 cases (only reached from m = ±3/2), and two for the other, as is also seen in the 
figure. Noteworthy is also the anti-crossing rule: even if some lines cross as B increases, no lines with the 
same rra-values ever cross. This is an effect of level repulsion between states of equal quantum mechanical 
symmetry (cf. Ref.Q pp. 292-295 and Fig. 10-2 therein.) 
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Figure 3: Definition of tlie polarization vectors. 
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Figure 4: The complete Zeeman pattern for the Hell hne at 4687 A (vacuum) consists of 146 components. It 
is here calculated for B — 2 Tesla at a viewing angle 9 = 40° relative to the magnetic field axis (cf. Fig.|3I). 
Only components with intensity greater than one (in the plot units) are included; the tt components are 
conventionally drawn with negative intensity. 
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